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1. INTRODUCTION 

A century-old history of calculation of asymptotics for solutions to Painleve 
equations (usually denoted Pj, j=l,2,...,6) as their variable x tends to infin- 
ity was started by pioneer works by Painleve, Gambier and Boutroux M. 
In 1980-1981 papers by Jimbo, Miwa and Flashka, Newell initiated the 
so-called Isomonodromy Deformation Method (IDM) which afterwards has 
given remarkable progress in this direction for P 1; P 2 , and P3. Successful 
attempts to apply the Whitham-like method connected with use of different- 
scaled variables were accomplished in papers ||. This way yields leading term 
of the asymptotic series in form of "modulated" elliptic function (Weierstrass 
p-function for Pi or Jacobi sine for P 2 ) such that parameters determining 
it depend on x in force of certain algebraic or differential equation. At last, 
the Whitham approach in its "classical" interpretation was applied to Pi by 
Novikov, Dubrovin, Moore, Krichever ||. These authors obtained certain 
ODE usually referred to as "Whitham "or "modulation" equations. In paper 
I]] solutions of Discrete Painleve- 1 eq. were investigated from this point of 
view. The main goal of this paper is to generalize ideas of the latter method 
and expand them to the remaining Painleve eqs. 

2. IDEOLOGY OF THE METHOD 
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Painleve eqs are known to be integrable nonlinear ODE's. Probably, 
the most explicit illustration of this fact is existence of commutative matrix 
representation 



d\Lj - d x Aj + [Lj, Aj] =o, j = 1, 2, .., 6, (1) 

where d x = Lj = Lj(u,u' ,x, A), Aj = Aj(u,u',x,X) are 2x2 matrices 
rationally depending on spectral parameter A so that j-th Painleve eq. u" — 
Pj(u,u',x) = is equivalent to flU). The matrices Lj, Aj are listed, for 
example, in 0. 

Introduce now variable X and change all variables x explicitly entering 
Lj, Aj to X : Lj = Lj(u,u',X,X), Aj = Aj(u,u', X, A). 

Proposition 1. Let e be some real positive number. Then system 

ed x L j -d x A j + [L j ,A j ]=0 (2) 
is equivalent to system 

d x X = e . . 

u" -Pj{u,u',X) = o [ ) 

Proposition 2. Solution of eq. (0) while e = and X = const is 

u (x) = fj(r + $;a*) (4) 

where r = xU ; fj are periodic functions which can be represented in terms 
of Weierstrass or Jacobi elliptic functions (see 0) for any of the six Painleve 
eqs. U = U(X), $ = $(X);~a > (X) are parameters determining the elliptic 
function fj. 

Suppose now e to be small and positive. We shall look for solutions of 
(|3]) as formal series in e : 

u{x) = u (x) + eu\{x) + ... (5) 

so that parameters determining the elliptic functions uq = fj obey some 
special ODE called "Whitham" or "modulation" system. Thus, we look for 
the leading term of (|) in form u (r, X) = fj (™ + HX); lt(X)) , where 
dxS = U. The Whitham system can be easily derived from eq. (0) by simple 
change d x — > Ud T + edx and averaging in r : 
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d x det Aj = a 22 d x ln + a u d x l22 - a^dxhi - a 21 d x h 2 (6) 

where Aj = (a mn ), Lj = (l mn ), m,n = 1,2; bar means averaging over period 
of the elliptic function 

Proposition 3. There exists only unique coefficient in the polynomial detA, 
whose dynamics in X is non-trivial. Denote this coefficient Fj. Therefore the 
Whitham system can be written out as a single first-order ODE on Fj. 

Corollary. The simplest way to find elliptic ansatz is to solve equations 

Fj = consti t X = const2 (7) 
Proposition 4. The Whitham system induces formal condition 

Uq - Pj(u ,u' 0} x) = 0(e)Q 

Obviously, the scaling transformation x — > ex induces d x X = e —> d x X = 
1 in (Q), which provides the following proposition: 

Proposition 5. The function u specified by (f|) and (||) yields leading term 
of asymptotic series for solution of appropriate Painleve equation as x tends 
to infinity. 

3. CONCRETE CALCULATIONS. 

Pi : u" - 3u 2 - x = 

det A x = 16A 3 + 4XA - F u where F 1 = (u') 2 - 2u 3 - 2Xu; d x det Aj = 
2u + 4A 

dxFi = —2u = —2r]/to = 2e\ + 2(e 3 — e-i)^, where K, E are complete 
elliptic integrals: 

K = J dz E = J J^fdz, k 2 = ei, e 2 e 3 are 

roots of the polynomial R^{z) = 4z 3 — g%z — g$, g 2 = —X, g 3 = —\F\. 
u = 2p(x + $; 5-2, ga). 

P 2 : u" - 2u 3 - ux = 0[] 

det A 2 = 16A 4 + 8XA 2 + X 2 - 4F 2 , where F 2 = (u'f - u A - Xu 2 ; 

^^This condition makes rigorous sense only for smooth function u. However, the method 
describes singular solutions as well. 

2 Case with additional constant can be investigated in the same way 
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d x det A 2 = 8X 2 + 2X+Au 2 => d x F 2 = u 2 = 2^-e m = X (± - ) , 

where k 2 = £r(X T VX 2 - 4F 2 ) - 1; 

m = \j-\( x =F V^ 2 - 4F 2 )sn C^yf-x ± V^ 2 - 4^2 + fc) , sn is Ja- 
cobi sine. 



P 3 : u" + \u' + sin w = 

A 4 det A 3 = ^ + ^fiF 3 + 1, where F 3 = - cos u; 

d x F 3 = -|(F 3 + cosU) = ^ (fc 2 - 1 + 2#) , A; 2 = (-F 3 ± sjFl^lj ' ; 

u = -2% log (fcsn (§ + 

P 4 : u " = ig! + 3|i + 4a;n 2 + u(a - 4d + 2a; 2 ) - where a, d are 
arbitrary parameters, det A 4 = — ^ — XA 4 — (d — a + X 2 )\ 2 — F 4 — p-, 
F4 = - i - Xu 2 - u(x 2 + § - 2d) + X(2d -a) -J, 
<9 x det A 4 = -A 4 -2XA 2 +w 2 +2Xw+a-2d <9 X F 4 = -w 2 -2Xw-a+2d 
and m = ±C/[C(C/s + $1) - ((Ux + $ 2 )] ± ^C($2 - $1) - X parameters 
U, $2 — $1 can be derived from (|7|). The Whitham eq. looks too large to be 
written out explicitly. 



p s : u" = (u'f - £ + & - fi 2 ^, n is a parameter; 



2uu'(fi - 4) + w 2 /i(/i - 8) 



2it(u-l) V 1 */ x 1 A" a A* 2(u-l) ' 

det A 5 = -4X 2 -^, F 5 = [K) 2 

d x det A 5 = -8X + jfy, where z = =s> d x F, = -Az; 

u = ^jf [p(Ux + $0 + p(Ux + $ 2 )] 

+ F Zi-ut) i((Ux + $0 - ((Ux + $ 2 )] + V, 
parameters [/, $1 — $ 2 , V can be derived from ([]). 



Pr- 



os 



u" 



1 x) 
z 2 (:r-l) 2 



1(1+1+^) ( M f _ (I + 1 + _0 u > 

2 \u u—x u—x J y ' \% x—l u—x ) 

-} 

(u—x) 2 _ 



a, P, 7, 5 are parameters; 
A [(A - 1) (A - a;)] 2 det A 6 = hk 2 \ 3 + F 6 X 2 + S t \ + S ; 
One can determine F$ and the elliptic ansatz from the following con- 
straint: 
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(u'f x 2 (x - l) 2 + u'a + (h - k 2 fu 4 + a 3 u 3 + a 2 u 2 + a x u + Q\x 2 = 0, 
where: 

§a = M 2 x(a;-1)(A; 1 + A; 2 )+Ma; [:r 2 (0 o + 6{) + :r(0 2 - 0i) - 9 - 9 2 ]-9 x 3 + 
9 x 2 , 

±a 3 = ihih - k 2 )(x + 1) + (3k 2 - k^xiOo + X ) + 9 + 9 2 ] - 2F 6 , 
a 2 = x 2 [(0 o + 00(00 + 9 1 + 4fci - 4A; 2 ) + Ak 2 (2k 2 - A*)] + 2x\l$k 2 - 6k ± k 2 - 
2k 2 - 

+0 2 0x + 2F 6 ] + Ak 2 (2k 2 - 2h) + (0 O + 2 )(0 O + 9 2 + Ah - Ak 2 ) + 4F 6 , 
iai =x 2 [(0o+0i)(2A; 2 -2A; 1 -0 o )+2A; 2 (A; 1 -2A; 2 )]+a:[(0o+^)(2A; 2 -2A; 1 -0 o ) + 
+2fc 2 (fc 1 - 2fc 2 )] - 2F 6 ; 

fci + A; 2 = -(0 O + 0i + 2 ), fci - A; 2 = 0^, a = |(0oo - l) 2 , (3 = -±0 2 , 
7 = |0 2 , 5=|(l + 2 ); 

d x F e = -2k 1 k 2 +z 2 (k 1 -k 2 )+9 2 (z 2 -k 2 )+ { -^^ - {zi+ ^- x) , where 
z = -^{u(u-l)(u-x)z 2 +[9 1 (u -x) + x9 2 (u - 1) - 2k 2 {u -l)(u- x)} z+ 
+k 2 (u -x-1)- fc 2 (0x + X0 2 )}, 
*i = -j^fcMu - 1)(« - x)z 2 + 

[(0i + 0oc)(w - x) + ^(w - 1) - 2£; 2 (w - l)(u - x)] z+k%(u-x)-k 2 (9 1 + 
x9 2 ) - k ± k 2 }, 

[0i(u - x) + x(0 2 + 0^)^ - 1) - 2A; 2 (w - l)(u - x)} z+k 2 (u-l)-k 2 (9 1 + 
x0 2 ) - xkik 2 }, 

^ u'x(x—l)—6o(u—l)(u—x)—6iu(u~x)—02u(u—l) 

2u(u—l)(u—x) 
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